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1. COEFFICIENT OF FRICTION 

Statement of the Problem. The friction Between metal surfaces 
completely separated "by a lubricant layer, or fluid friction, is the 
reaction of the layer to the motion of the surface. Experiment 
and theory show that fluid friction does not depend, or depends only 
to a slight extent, on the magnitude of the load on the surface and 
is determined to a very large extent by the viscosity of the lubri- 
cant and the velocity of the relative motion. 

The usual method of expressing the friction force as a function 
of the load through the coefficient of friction is for fluid fric- 
tion quite an artificial one and its adoption in the lubrication 
theory is merely a concession to the concepts of friction based 
on the laws of Coulomb. 

The numerous experimental and theoretical investigations of 
the coefficient of friction of a Journal bearing all express the 
coefficient as a function of the parameter T]m/k, where rj is 
the absolute viscosity of the oil, to the angular velocity of 
the shaft and k the pressure. This magnitude which plays a 
large part in the discussion that follows we shall denote as the 
operating parameter of the bearing. For practical computations it 
is more convenient to make use of the related magnitude T]n/k, 
where tj is the speed of the shaft in revolutions per minute, k 

the pressure in kg/cm 2 on the projection of the bearing. For 
briefness Tjn/k will be denoted by the symbol X. The latter is 
connected with the magnitude ijco/k by the following relation 

\ cen, fclpolses rpm _ g^6Q • iQ ^jjjg, sec/m 2 cu l/sec 
kg/cm 2 k kg/m 2 

Petroff in 1883 proposed for the coefficient of friction of 
parallel surfaces separated by a lubricant layer of height h and 
moving relative to each other with velocity v the following expres- 
sion: 
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f = L . 31 
h k 

which follows immediately from the egression of Newton for the 
force required for the displacement of a viscous fluid between two 
parallel surfaces of solid bodies. According to the law of Newton 
the shear force T of a viscous fluid is proportional to the area, 
the viscosity of the fluid, and the velocity gradient across the 
fluid layer: 


T = lb t) 


dv 

dy 


( 1 ) 


where 1 is the length and b the width of the surface. 

The shear force thus depends on the velocity profile across 
the oil layer. If the oil in the clearance is carried along only 
by the adhesion to the moving surface and the viscosity force the 
velocity gradient is constant and equal to v/h, the velocity 
profile is linear and equation (l) assumes the following form: 

T - lb ij ?- (la) 

h 

The coefficient of friction or the ratio of the frictional 
force T to the load on the bearing surfaces P = klb In this 
case is equal to: 


f . i m 
h k 

For a cylindrical shaft rotating concentrically with the 
angular velocity o> in a bearing with diametral clearance A the 
friction force by the law of Newton is equal to: 

T = it dl ^ 

A 


and the coefficient of friction 

f = 1 s 1 (2) 

kid A k \[/ k 

where \j / = A/d Is the relative clearance of the bearing. This is 

the formula of Petroff for a Journal bearing. It was derived for a 
clearance of uniform width and takes Into account only the friction 
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of the viscous displacement of the oil and is therefore applicable 
only to the determination of the friction of lightly loaded and 
concentrically rotating shafts. ... . ... 

If a load is applied to the surfaces then, in addition to the 
reaction of the viscous displacement, the surfaces are subject to 
still another reaction of the oil flowing out under the pressure in 
the direction of the motion (fig. 2). For cylindrical bearings 
the zone of pressure constitutes generally 90 to 120°; In this zone 
forces are developed which increase the friction of the shaft above 
the value determined by the Petroff formula. 

In determining the simultaneous reaction of the viscous dis- 
placement and potential flow of the oil in the loaded zone, assumed 
as 120°, Gurnbel (reference 2) expressed the coefficient of friction 
as a function of the relative eccentricity of the shaft in the 
following manner: 


f = 1.7 Vl^X (3) 

where X is the relative eccentricity equal to the ratio of the 
absolute eccentricity of the shaft e to the radial clearance 
8 a A/2 (fig. 3). The relative eccentricity, as follows from the 
equation of discharge of Reynolds is a function of the magnitude 
T|oj/k\J/2 that is 


<’’> -gj? M 

Substituting in equation (3) for 1 - X the algebraic expression, 
represented in figure of the relation between X and tjcjo/ k^, 
obtained by the numerical integration of a number of typical cases 
Gumbel obtained for the coefficient of friction the following expres- 
sion: 


f . 1.7c tfmT 

-Hie coefficient e that takes into account the Inaccuracy of the 
algebraic expression for the relation between X and tjcjo /Hk 2 Is 
shown in figure 1 as a function of T)cc/k\l? . 

The eccentricity of the shaft in the bearing depends not only 
on the magnitude rjco/ky 2 but also on the length of the bearing. 

The smaller the ratio l/d the more easily doeB the oil flow out . 

Of the clearance and the deeper must the shaft be set in the bearing 
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In order to sustain the given load. To take into account the finite 
length of the ■bearing it was necessary for Gumbel to introduce in 
the above derived expression the coefficient of length of the "bear- 
ing obtained from the working up of the tests of Lasche (reference 3) . 
The final expression for the coefficient of friction according to 
Gumbel has the following form: 

-c 



where c* = yk d/l + 1 is the correction for the finite length of 
the bearing. Assuming the value of e equal on the average to 1 
and choosing the ratio l/d = 1 Falz (reference 4) obtained the 
following widely known formula: 

f » 3.8^ (6) 

The above formula gives satisfactory results for partial bearings 
(without upper cover) with a bearing arc of about 120° but is unsuit- 
able for full Journal bearings particularly those with forced feed 
since the formula doeB not take into account the additional fric- 
tion of the viscous displacement of the oil over the unloaded two 
thirds of the bearing. The formula may be used for the approximate 
estimate of the friction in full Journal bearings only for small 
values of Tp/k where the additional, friction of the viscouB 
displacement is relatively small and the reaction of the flow from 
the loaded zone is the predominant factor. Moreover the range of 
values of rjcu/k in which formula (6) gives satisfactory results is 
so near the critical value of Tp/k corresponding to the rupture 
of the oil film that it is practically Impossible to use 
formula (6) for full journal bearings. 

Test Data on the Coefficient of Friction. Test on the coef- 
ficient of friction generally determine the latter as an exponential 
function of tj, cd, k. Ulmer (reference 5 ) on carefully working up 
almost the entire literature on the coefficient of friction 
represents the latter in the following form: 


f = K &L 
hP 

where the exponents r, m, n, vary within the range 0.3 to 1 
and K is a constant that takes into account the size of the clear- 
ance and other design factors. The rather complicated method proposed 
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by Ulmer for determining the coefficient of friction is har dly 
applicable on account of its complexity and absence of any 
rational basis. We present the results of the experimental investi- 
gations of the various authors considered by Ulmer (table l) in 
liable 1 . 

Fomin (reference 6 ) in conducting tests with the crankshaft of 
an airpl an e engine found that r varies in the mean between the 
limits 0.3 to 0.4, m = 0.5 to 1 , n = 0.6 to 1 . 

McKee (reference 7) investigating the coefficient of friction 
and its dependence on the clearance and length of bearing expressed 
the results of his tests by the following equation 


f 


2L. IE. 

\J 1 k 


+ 0 a 


the first term of which represents the formula of Petroff, a is a 
constant equal to 0 . 002 , 0 is a correction coefficient depending 

on the length of the bearing. 

Coefficient of Friction of a Full Journal Bearing. The velocity 
distribution in the oil layer of a bearing working in the region of 
fluid lubrication with a clearance completely filled by the oil is 
schematically represented in figure 2 together with the diagram of 
the pressure over the circumference of the bearing. Over the entire 
circumference of the bearing the oil is carried along by the moving 
surface of the shaft and the velocity diagram across the oil 
layer had the form of a triangle the maximum ordinate of which iB 
equal to the peripheral velocity of the shaft. The resistance of the 
oil, carried along by the motion of the shaft, per unit area is 


T 




(7) 


In the narrow part of the clearance due to the incompressibility 
of the oil a region of pressure is formed from which the oil flows 
out in the direction of motion and toward the ends of the bearing. 

The latter type of flow has no effect on the friction since the 
frictional force is directed at right angles to the peripheral 
velocity. The flow in the direction of motion however (along the 
X axis) gives additional friction. To the triangular velocity 
profile there is here added the parabolic profile of the potential 
flow, the equation for which, for a system of coordinates placed 
at the center of the oil layer, has the following form: 


6 


naca m no. 1165 


v 



dx 


where dp/dx is the pressure gradient along the X axis. 


(8) 


The first derivative of the velocity with respect to y is 
equal to 

dv y dp 
3y “ T) 3x 


The friction force per unit area for the boundary surfaces, that is, 
for y = + h/2 is according to the law of Newton equal to 


dv li dp 
11 dy = 2 dx 


(8a) 


The total friction of the hearing is made up of two components: 
of the viscous resistance of the oil T* carried along by the motion 
of the shaft over its entire periphery and of the reaction T" of 
the flow of the oil from the loaded zone along the direction of 
the X axis: 


T = T* + T" * dl (jit + aT 4 '), 

where a is the loading angle of the* bearing. The total coef- 
ficient of friction is the sum of the individual components 

f . + - f + f 

p p 

T* and f 1 . We shall determine the value of the first form of 
the friction according to equation (la) . We transform the linear 
coordinates of this equation into polar for which purpose we make 
use of figure 3 showing the shaft eccentrically placed in the bearing. 
We take the origin of coordinates on the line 0-0 connecting the 
centers of the shaft and bearing. For any angle cp the height 
of the clearance h can then be represented by the following 
egression: 


h = r \j/ (1 + X cos cp). 


(9) 


where 

r is the radius of the shaft 

(jz the relative clearance, equal to A/d 

X the relative eccentricity, equal to 2e/A 
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The total friction force over the entire periphery of the shaft is 
obtained by integrating expression (10) between the limits 2* to 0: 

T » „ irjyp 2 n d{ P r -S 

\i/ J 1 + x cos <p ^ yi-x 2 

0 


The coefficient of friction f* is equal to 

f* . _I" =JL3fiJ —l »XM (X) * (11) 

kid ^ k yl - I 2 \J/ k 

VTe have thus obtained the formula of Petr off corrected by the 
factor —J. .. that takes account of the eccentricity of the shaft 
Vl - X 2 

in the bearing. For X = 0, that is, for a central setting of the 
shaft, formula (11 ) becomes the formula of Petroff . 

Figure 4 for the function (X)* shows that the correction term 
differs to an appreciable degree from 1 only for values X>0.5 
and that therefore the formula of Petroff is sufficiently accurate 
over the large range of practically occurring small eccentricities 
corresponding to high values of the ' parameter qco/k. 

For convenience in using formula (ll) we replace the eccentri- 
city X by the magnitude T)<n/]fy 2 of which it is a function. The 
relation between these magnitudes is given by the flow equation of 
Reynolds which in the interpretation of Gumbel (reference 9 ) for 
the cylindrical bearing for a mean value of the loading arc of 120° 
leads to the relation shown by the curves in figure 5- Tbe lower 
curve refers to a bearing of infinite length. For bearings of 
finite length Gumbel introduced in the expression of the relation 
between X and r)u>/k\)F the correction factor c = 1 + d/l 
obtained by working up the test data of Lasche and his own data: 

The author found that the function (X) for a bearing of infinite 
length, as shown in figure 4 is satisfactorily expressed in terms 
of rjaykty 2 by the following empirical equation: 

/— \ 4 , ' « . 0«X /ia> 
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For a bearing of finite length it is necessary in equation (12) to 
introduce the correction factor c. 


Substituting equation (12) in expression (ll) we obtain 

f (1 + ] -f ^-+ o.314c 

\ 1 / k \ ^ / \1 / k 

\ / 


(13) 


that is, the first partial coefficient of friction is equal to the 
friction coefficient of Petroff for the concentrically rotating 
shaft plus a certain magnitude, constant for a given bearing, that 
expresses the effect of the eccentricity of the shaft in the bearing. 
It is very important to notice that the finite length of the 
bearing is expressed only through the constant term of equation (13) . 
The smaller the value of l/d the greater for a given value of 
qoj/k the eccentricity of the Bhaft and the greater the coefficient 
of friction. 


T" and d" . We found previously (8a) that the friction force 
per unit area due to the flow in the loaded zone of the bearing is 
equal to: 


T" = 


h dp 
2 dx 


or in polar coordinates 


t" (1 + X cos cp ) dp 
~ v 2 d<p 


(HO 


The equation of Reynolds, transformed for a cylindrical bearing, 
glvea for the pressure gradient over the circumference of the bear- 
ing dp/dcp the following expression: 

dp r oi (l + X cos <p) - (1 + X cos cPn) 
d.q) = &TI ^"2 (1 + X cos <p)3 


where cfo is the angle corresponding to the point of minimum clear- 
ance. It is simpler however to obtain the pressure gradient, with 
an accuracy sufficient for practical computations, if we make use 
of the test data on the limits of the bearing film and the % character 
of the pressure change over the loading arc. 


Numerous experimental investigations all give the following 
picture of the pressure distribution in the plane of symmetry of 
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the full cylindrical hearing under constant load. The pressure, 
equal to zero over the entire periphery of the hearing (if the oil 
is fed without pressure), begins to increase over 90 to 1200 
from the line connecting the centers of the shaft and hearing, and 
increases with an almost constant gradient up to a maximum value 
lying near the point of narrowest clearance and constituting on 
the average 2.5 to 3 k. Thereafter the pressure drops sharply to 
zero or even somewhat helow. A typical pressure curve measured 
hy Bradford and Grander (reference 10) is shown in figure 6. 

Without great error this picture may he replaced hy the scheme 
shown in figure 7. The pressure here increases linearly over an 
arc of 90°, attains the value 2.5 k at the point of minimum clearance, 
after which it drops very steeply. The assumed scheme is suffici- 
ently near the true conditions, the approximation being closer the 
smaller the value of rjco/k and, as we shall see helow, our results 
are of importance precisely for these values of qcc/k. This 
assumption, very much simplifies the computation of the coefficient 
of friction hy equation (14) . .The pressure gradient dp/dp 
becomes constant and equal to Pttwt / 8 - = sk/a where a is the 

loading arc which in what follows we shall assume equal to 
jt/2, s = 2.5 - 3.* For every other equivalent section of the hear- 
ing not lying in the plane of symmetry the magnitude dp/d<p is 
equal to 2p m /« where is the pressure at the point of 

minimum clearance in that plane. 

The pressure distribution over the length of the hearing (the 
Z axis) according to test data agreeing with the theoretical con- 
siderations (reference 11) is represented in any meridional section hy 
a parabolic type curve with almost constant exponent m varying 
between 2.2 to 2.7. The maxi mum of the curve is on the axis of 
symmetry of the hearing as shown in figure 8 giving the pressure meas- 
urements of Nucker (reference 12) over the length of the hearing. 

The pressure p may he expressed as a parabolic type function 
of the maximum pressure = sk: 



"The computations made hy the author for a = 9° to 120° show 
that the results given helow are valid for the entire practically 
occurring range of values of a. 
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and, the pressure gradient at any point of the loaded zone of the 
"bearing may "be expressed as follows: 

dp ^ 2pm = 2sk 
dcp it it 

The different points of the surface along the length of the 
bearing thus experience a different reaction t_" as a function of 
their distance from the plane of symmetry of the bearing. 



On. the loaded surface of the shaft consider an element of width 
rdcp and length dz (figure 9). The force dT" acting on this 
element is equal to 


dT" « t" rdcp dz = rdcp dzV - & - + ■ x , cog ft) £2 

2 dtp 

= rdp dz* Sk± L.. g -°!L±1 

2 n 


-a 


(15) 


To determine the total force T" acting on the shaft we integrate 
equation (15) twice, once over the length of the bearing, that is 
between the limits ±L/2. and a second time over the length of the 
bearing film, that is between the limits it to it/2; 


T" = 


1 . 

it + 2 




X cos <p) 


2 2 


1+1 



m 


d p dz 


= _m_ {1 _ 0#6 ^ J)# 

2 m + 1 


The friction coefficient f" is equal to: 


f" - = p . ..- m - (1 - 0.64 X)V s -2_ (x) 

kid 4 m + 1 4 m + 1 L 


Assuming s = 2.5 and m = 2.5 we obtain 

f" = 0.4^(X) A 
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The function (X)^ expressed In terms of t\(d/TpI^ has the follow- 
ing form 



1 




1 + 2.5 



The Relative Magnitudes of f * and f w . The weight of each of 
the component parts of the total coefficient of friction is shown 
in tables 2 to 4. Those taDJ.es give the values of f* and f" within 
the range = 0.01 to 100 which for the relative clearances 

\l/ = 0.001, 0.002 and 0.0005 corresponds to the values of the 

operating parameter: 


1 — 10 . 100 000 centipolses X rpm/kg/cm 2 


The double line in the tables divides off those values of Tp/kxJ? 
which in all probability axe outside the brakedovn limits of the 
oil film. Examination of tables 2 to k shows that the term f" 
affects the total value of the coefficient of friction only in the 
range of small values of rjo/kl / 2 which probably lie outside the 
limits of fluid friction. For large values of rjo/k^ 2 the effect 
of the factor f" iB not large. Computation shows that f" may 
without great error be assumed constant and equal to its mean value 
in the most important interval of Tjoi/tol/ 2 : 


f" = 0.25\|r 


The error arising from this assumption in the least favorable case 
does not exceed 5 percent.* 

Total Coefficient of Friction. Under the above made assumption 
the total coefficient of friction takes the following form: 


f - f » + f”^- 1 ^ + 0 .31W + 0.25^ - f ^ + 0.55^ ( 16 ) 

y k * 

The ratio l/d, as follows from the previous considerations, 
affects only the constant term of equation (l6) as is fully 
confirmed by the above mentioned tests of McKee . According to 


*The relative weight of f" in the general value of f has 
been computed on the assumption, of the same viscosity of the oil 
over the loaded and nonloaded bearing. Actually the temperature 
of the oil in the loaded some is always somewhat higher than in the 
remaining part due to the additional friction and the relative value 
of f" is somewhat smaller than that given in tables 2 to h. 
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these tests for l/d> 1 the consteuit term of the expression for the 
friction coefficient must "be corrected hy the factor a which may 
he expressed hy the following equation: 


a 



1.5 


For values of l/d> 1 the correction factor o is approximately 
equal to 1. 

With the correction for the length of the hearing the friction 
coefficient takes the following form 

f 'f* + °-55(x) 1 ' 5 * (17) 

or in practically suitable units: 

f = 3.36 • 10-9 d x + 0.55^) 1,5 A (18) 


where X = rjn/k 

d the diameter of the shaft in ram 
1 the length of the hearing in mm 
A the diametral clearance in ram 
q the. viscosity of the lubricant in centipoises 

n the rpm 

k the projected unit hearing load in kg/cm 2 

The correction (d/l)^-*5 is introduced only for values of l/d less 
than unity. 

Comparison with Test Data. Figure 10, to logarithmic scale, shows 
the curves of the coefficient of friction as systematised hy A. Wewerka 
(references 16) according to the test results of a number of authors 
(references 13, l 1 *-, 15, 17) the formulas of Petroff, Gumbel and 
formula (17) of the present paper. The curves are plotted in the 
coordinates proposed hy Howarth, f /J / being laid off on the ordinate 
axis against ticd/Jo^ on the abscissa. 
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it 

The curves on figure 10 show clearly that the formula of Gumbel 
satisfactorily agrees with the test results for small values 
of T^cu/k^ hut gives a large error for high values of this factor. 
The formula of Petr off, on the contrary, agrees very well with 
test results for large values of Tj<e/k\|£ as was to he expected 
from the structure of the formula, hut gives too small, values far 
f for small values of qo>/kv|^ far which the eccentricity of the 
shaft and the friction in the loaded zone have a large effort. 
Formula (17) satisfactorily agrees with test results over the entire 
range of values of rjco/k^ practically encountered. 


The dotted straight line drawn through the mean values of the 
coefficient of friction may serve for an orientating estimate of the 
coefficient of friction in the range r\(a/'k^ = 0.5 to 100. The 
equation for this line is the following: 


f * 5 fm\ °* 8 5 

* 0.7 V k j 


(19) 


or in practical units 

f - 12 • 10- 8 ^°-T X °- 8 5 ( 19 ) 


where X = rjn/k, 

tj the viscosity in centipoises 

n the rpm 

k the projected unit hearing load in kg/cm 2 

d the di ame ter of the shaft in mm 

A the diametral clearance in mm 

General Character of Formula (17 ) . We may note that formula (17)# 
derived on the assumption of constant load on the hearing, may he 
extended with a certain degree of reliability to other loading cases. 
Formula (17) without any reservations is applicable to the case of a 
shaft loaded by the rotating vector of a centrifugal load. This case 
represents the rotating system of a hearing under constant load with 
the characteristic feature that the load k is a function of the 
square of the speed and the coefficient of friction is a function of 
the speed and viscosity: 


Ik 
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We may note that in thiB case the operating parameter drops no 

with increased speed, in contrast to the case of constant load, and VL 

the coefficient of friction decreases with increased speed. For an 
impact load and for a load that varies in magnitude and direction 
formula (17) may serve as an approximation which is the more accurate 
the larger the value of qa)/k\^. The character and magnitude of the 
friction in the unloaded zone are here the same as in the case of 
constant load hut the friction in the loaded zone has a different 
character . The curve of pressure in the loaded zone in all proha - 
hility approximates a parabola the axis of symmetry of which for small 
speed of rotation passes through the line connecting the centers of 
the shaft and hearing and at large speed is displaced against the 
direction of motion. The oil in the loaded zone flows not only 
against the direction of motion as in the case of the constant load 
hut also along the direction as a result of which a considerable 
decrease and even the total disappearance of the reaction of this 
flow on the shaft may he expected. The probable diagram of the 
velocities along the periphery of the shaft is shown in figure 11. 

On the other hand the relation between the load and the eccentricity 
of the hearing for an impact load and a load varying in magnitude 
and direction is other than for the case of a constant load and is 
determined not only by the parameter T^oj/k hut also by the rate of 
increase of the load and the rate of change of its direction. The 
eccentricity is in this case over a large part greater than in the 
case of a constant load. 

Thus the variable and impact loads while leaving unchanged the 
first term of equation (17) change the constant term, the effect of 
the latter being in opposite directions. On the one hand the constant 
term decreases as a result of the flow in each direction in the 
loaded zone and on the other hand it increases as a result of the 
eccentricity. A detailed a n alysis of the friction for an impact and 
variable load constitutes a separate independent problem which the 
author hopes to consider another time. 

Formula (17) may serve for an approximate estimate of the coeffi- 
cient of friction in the case of an impact and variable load not only 
because of the opposite effect noted above the friction on the con- 
stant term of this equation but mainly because of the small relative 
value of this term in the practical interval of values of rjco/k . 

As shown by tables 2 to k the constant term of equation (17) even for 
values of Tjm/k near the critical does not exceed 25 percent of the 
total coefficient of friction. Actually bearings always work with a 
certain factor of safety and the magnitude of the constant term is 
generally very small in comparison with the first term. 
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APPLICATION CF THE OBTAINED RESULTS 


1. kv as a Measure of the Thermal S trees of the Bearing. 

Formula (17) answers clearly the question often raised as to the 
applicability of the magnitude kv as a measure of the heat 
generated in the hearing. This method of estimating the thermal 
stress implicitly assumes the constancy of the coefficient of friction 
for all operating conditions and for various structural design factors 
of the hearings and at first view has no basis in fact if the friction 
theory based on Coulomb is rejected. Formula (17 ) shows that this 
is not at all the case and that under certain conditions the heat 
generation is actually proportional to kv. Nevertheless a careful 
analysis of the problem makes it necessary finally to reject the 
magnitude kv as well as any other magnitude even for a rough 
estimate of the thermal stress of the bearing. 


The heat given out in unit time over unit surface of the 
bearing is expressed as a function of the coefficient of friction 
in the following manner: 


k cal 
sec • 


Pvf kv / « t|cd 
I d 427 “427 \k k 


1 

wr 


2n 


TJV C 


0.55 <niy 


"i A 

( 20 ) 

/ 

of two parts: 

One 


part proportional to tjv 2 and the other to kv. The relative values 
of the two parts are shown in tables 5 and 6 where the magnitudes 
of the component parts of equation (20) are given as a function 
of ijo>/k\l£ far l/d = 1 an&\|/ « 0.001 and 0.002. 


Tables 5 and 6 show that the magnitude kv predominantly 
determines the general value of the work of friction for small 
values of Tpi/ki^ lying probably beyond the limits of fluid fric- 
tion with the exception of those cases where the bearing has a 
small clearance and l/d is considerably less than 1. In the 
range of fluid friction the general value of R is greatly 
affected by the term proportional to tjv 2 and for large values 
of rjcu/kV 2 this term completely determines the value of R. 
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Thus the application of kv as a measure of the heat generation 
and as a basis for a comparative estimate of the thermal stress of 
the bearings may be justified for small values of Tp/kty 2 approaching 
the critical value. A second condition is the equality of the 
relative clearance i|/ and the ratio l/d for the bearings com- 
pared. For most bearings •working with a sufficient safety factor 
the use of kv gives a large error. This however is not the main 
point. The determination of the load -carrying capacity of the bearing 
from the value of the heat generated involves a methodological 
error. The load-carrying capacity of the bearing is determined 
for a given load, speed and kind of oil by the temperature of 
the bearing which inturn is determined by the mutual interaction 
on the one hand of the heat generated and on the other the heat 
dissipated. The latter depends on a whole series of factors, for 
example, the size of the clearance, the feed pressure of the oil, the 
cooling surface, the number and arrangement of the oil grooves and 
the operating conditions characterized by the parameter rfu/k. 

On the basis of the above considerations the magnitude kv or 
any other magnitude must altogether be dispensed with for even a 
comparative estimate of the thermal stress of the bearing and the 
only rational way is to determine by computation the heat balance of 
the hearing and the effect of the heat generation and heat desslpation 
on the load -carrying capacity of the bearing. This will be done below. 

2. Floating Bushes. Expression (17) sufficiently explains the 
advantages, often confirmed in practice under certain conditions, 
of floating bushes thus throwing doubt in the validity of the theories 
of the constancy of the coefficient of friction as shown below. 

In a stationary bearing the heat generated is R = const Pvf . 

Now consider a floating bush. If the coefficient of friction is 
considered constant it is assumed for simplicity that the 
floatin g bush rotates at a speed equal to half that of the shaft 
then on the two sides of the bush the amount of heat generated is 

E' = 2 const Pf v/2 

that is, the floating bush gives no gain In the heat generated. We 
shall now apply equation (17) assuming that the bearing works in 
the region of high values of Tpo/k so that the constant term in 
equation (17) may be neglected. The heat then generated in the 
stationary bearing is equal to 

B’ = const qv 2 

The heat generated on both sides in the floating bush (assuming the 
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clearances the same) Is 

R* = 2 const T] v^/4 

that is, the floating hush lowers the heat generated by half. Two 
concentric floating bushes applied in certain especially high speed 
bearings lower the heat generation to one third. The smaller the 
value of T]to/i: and the greater the weight of the constant term of 
equation (17) the less advantageous is a floating bush. If it is 
also remembered that the floating bush decreases the speed of rota- 
tion and hence the load bearing capacity by two the range of favorable 
application of floating busheB becomes entirely clear, namely, the 
range of high values of Tjco/k of high speed, low load shafts working 
with a high factor of safety. In this case the floating bushes 
provide a technically rational means for lowering the temperature 
of the bearing though it is true at the expense of a lowering in the 
safety factor. We shall not here touch upon the other extreme of 
the favorable application of floating bushes, namely, in the range 
of half dry and half fluid friction. 

3. Cut-away of Bearing and Shaft. A second method that has 
long been used in practice for decreasing the friction consists in 
increasing the clearance in the nonloaded part of the bearing as 
shown in figure 12a, or using a partial bearing as shown in figure 12b. 
For centrifugal load the same result is obtained by a cut-away of 
the shaft as shown in figure 12c. The advantages of the methods of 
figure 12a and 12c 1 b twofold. In the first place the friction in the 
nonloaded region becomes negligibly small and In the second place 
the flow of oil is increased. This will be considered in detail in the 
section on the oil flow. 

For partial bearings and cut away shafts and bearings the friction 
may be computed by the formula of Gumbel if the loading arc does 
not differ too greatly from 120°. Figure 10 permits the immediate 
determination of gain from the cut away of the shaft and bearing. 

For this purpose it is sufficient to compare the straight line rep- 
resenting the equation of Gumbel -Falz with the curve represented by 
formula (17) for a definite value of qto/k\|/ 2 . The gain is greater 
the greater the value of qfo/k^. Thus according to figure 10 for 
rjoa/k \}/ 2 = 5 the heat generated is reduced by half in comparison with 
the normal bearing and for T)co/k\l£ = 20 it is decreased 3*5 times. 

For values T]o>/k4^<l the cut away of the shaft and bearing has no 
significance. 

2. FLOW OF OIL FROM THE BEARING 

The flow of oil from a full Journal bearing with froced feed 
consists ©f two parts: (l) the flow q* from the loaded part under 
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the action of the pressure developed In this zone, (2) the flow q" 
under the action of the forced feed. 


Determination of q*. For the solution of this problem we shall 
make use of the scheme, already applied by us, of the variation In 
pressure in the loaded part of the bearing, that is, we assume that 
the pressure in the loaded zone varies linearly over an angle of 90 °. 


Die velocity of the oil along the Z axis (fig. 13) in the 
clearance of a bearing of height h according to equation ( 7 ) will 
be 



The volume of oil flowing per unit time through the clearance of 
width b is equal to 


<L 



vdy 



dp bh^ dp 

dz ^ = 12 tj dz 


( 21 ) 


09 

ro 

V£> 


We shall transform equation (21) into polar coordinates. We again 
turn to figure 3 , taking on the surface of the shaft an element of 
width rd<p and height h - r (cos <p) . The flow through this element 
in the direction of the Z axis is equal to 

rd <p (1 + X cos <p)3 r^\i/ dp . . 

dq ■ — (22; 

* 12tj dz 

As before we assume that the pressure along the axis of the bearing, 
that is, along the Z axis varies parabolically with exponent m. 
The pressure gradient dp/dz at the end of the bearing is equal to 


dp mp 
dz ~ 1 
2 


(23) 


where p is the pressure in the plane 
The variation of p over the range of 
of our mathematical assumptions may be 
equation : 



of symmetry of the bearing, 
the bearing angle on the basis 
expressed by the following 



J 
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where jK mx « sk Is the pressure at the point of minimum clearance 
In the plane of symmetry of the hearing., 0 Is a coefficient fluc- 
tuating within the limits 2.5 to 3 and a the loading arc. Par 
a = */2 we have 


P “ 



(24) 


Substituting equation (24) In expression (23) and the latter in equa- 
tion (22) we obtain the following expression for the flow through the 
element of the clearance in the loaded part of the bearing: 


dq* 



6t]1 




1 + I cos <p/3 d<p 


(25) 


The total flow over the entire loading arc is obtained by integrating 
equation (25) within the limits n - tr/2: 

,• - fa - + I oos „>3 ^ 

orjl JL\ rt / 

2 

r 4 \|r 3 msk/ 0 , „ , . , , 2 _ i -,_3\_ r 4 \|f 5 msk m 

= — ^2.1 - 4.3 X + 3.7 Z - l.lj 6tj1 w o 


For the two sides of the bearing 

r 4 \lrmsk 

«’ ‘ it (I) ° (26 > 

The magnitude k is expressed as a function of the relative eccentri- 
city X by equation (3) : 


k = 


2EL 

* 2 (XV 


where (X)* is the function represented in figure 1. Equation (26) 
in this case takes the following form: 

1 „ (I)q 

q " 31 (I)* 


( 27 ) 
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The magnitude (X) Q / (x)* is shown plotted in figure 14 as a function 
of X. As shown hy figure 14 the flow from the loaded part of the 
hearing decreases with decreasing eccentricity of the shaft in the 
hearing and becomes equal to zero for a central position of the shaft. 
This was to he expected since For such position the hearing can not 
carry any load and there is no pressure in the oil film. 


For X « 1 the flow is infinitely great since according to the 
Reynolds equation the pressure in the narrowest part of the clearance 
becomes infinitely great. Such a position never occurs in practice 
since long before this half -dry and dry friction arises as a result 
of the contact of the surfaces of the shaft and hearing. 

The magnitude (X) Q /(X) # may with an accuracy sufficient for 
practical purposes he expressed as a function of hy the 

following empirical equation: 


(2lo 

(«* 


0.25 


1 tjw\0.6 

7 i z¥j 


(28) 


where c is the correction factor of Gumbel for length of hearing. 
Substituting expression (28) in equation (27) we obtain 


c^2«2 xnsco^l 4- 

is i ^°- 6 


Substituting 



rfn 

» — , s * 2.5 m =2 .5 
30' 


0 »6 

aru * replacing (i + d/l) by an approximate formula from the binomial 
expansion we obtain 


where 


q’ 


0.0034 



cf 



0.6 



(29) 
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The flow per rotation 1 b equal to 

qn* - 0.2 £ - ^ 2 * 2 gt (30) 

^>) 0.6 

Reducing equations (29) and (30) to practical units we obtain 


. ^ „ d A o*n 

q cm 0 /soc ■ 0.8 ^ 

^0«o 

(31) 

_ ..0 .8. 2 .2 _ 

qn* cm 3 /rotation = 50 

(32) 


x o.S 


where X ■ r^n/k 

d the diameter of the shaft in ran 
A the diametral clearance in mm 
T) the viscosity of the oil in centipoises 
n the rpm 

k the projected unit bearing load in kg/cm^ 

1 the length of the bearing in ran 

— (1 + 0.6 d/ 1 ) 

1 

Equations ( 29 ) to (32) show that the flow of oil in the bearing with- 
out forced feed to a large extent depends on the clearance, on the 
ratio 1 /d, is almost directly proportional to the diameter of the 
shaft and inversely proportional to the parameter X. For a given 
bearing the oil flow is a function of X only and follows the law 
represented in figure 14. 

The Bering as a Stable System. Figure 14 illustrates the notable 
property of a slider bearing observable not only for a bearing without 
oil circulation but as we shall see later also for a bearing with 
forced circulation. This property is namely the fact that the bearing 
working in the range of fluid lubrication constitutes a system in 
equilibrium tending to maintain a definite value of the parameter qco/k 
against all disturbing factors. 
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Let us compare figure 14 with the typical curve of the friction ^ 

coefficient plotted against X (fig. 15) • We shall assume that the 
hearing works under a steady condition for example at X = 1500. 

We shall assume further that some disturbing factor enters in the 
operation of the hearing, for example, an increase in the tempera- 
ture due to local friction. In this case the viscosity of the oil 
drops and X decreases. This produces on the one hand a decrease 
In the friction according to figure 15 and on the other hand an 
increase In the flow according to figure 14. The simultaneous action 
of these factors lowers the temperature and reestablishes the pre- 
vious value of X. An analogous picture is obtained for a decrease 
in X due to and increase in k. 

The reverse phenominon occurs on increasing X. Let us assume 
for e xam ple that X increases due to an increase in the speed. In 
this case the friction according to figure 15 increases, the flow 
simultaneously decreases and the temperature of the bearing rises 
as a result of which the viscosity of the oil decreases and the X 
tends to drop. 

The automatic character of this process is due to the property a 
of the oil of changing its viscosity with the temperature and the 
unique dependence of the flow on X. In this peculiarity may be 
found the secret of the notable stability and noncapriciousness of 
behavior of correctly designed and constructed bearings. The essential 
condition for this automatic stability is that the bearing works with 
a sufficient safety factor and the time fluctuation of X should 
not carry it beyond the critical values at which a breakdown of the 
oil film occurs. 

Range of Application of Formulas (29) to (32) « Equation (29) 
to (327 are applicable to the majority of technical cases of the 
bearings without forced feed lubrication or for a feed pressure not 
greater than 0.1 to 0.5 at. The essential condition for the applica- 
bility of formulas (29) to (30) is the absence of oil groves in the 
loaded zone and the rigidity of the bearing and shaft. Such grooves 
even though they do not extend to the ends of the bearing may increase 
the flow many times. However modern bearings only rarely have grooves 
in the loaded zone. 

The deformation of the bearing or shaft under the action of the 
load may change the law of pressure distribution along the axis of 
the bearing and the quantity of oil flowing out. 

Notwithstanding the fact that formulas (29) to (32) were derived 
on the assumption of constant load they may without reservations be 
applied to the case of a rotating centrifugal load which may be trans- 
formed into the scheme of constant load. It should be noted that the 
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unit load, k la here a function of the square of the speed, and. on 
increasing the speed, the parameter X, in contrast to the case of 
constant load., decreases and therefore the flov Increases. This 
has long heen observed on systems in which the centrifugal load 
predominated, for example, on crankshafts of Internal combustion 
engines, etc. 

In the case of an Impact or variable load equations ( 29 ) to (32) 
may serve only as a first approximation since the law of pressure 
distribution along the loading arc in this case differs from that 
assumed as the basis of these equations which as a result give de- 
creased values of q ' . 

Oil glow for a Bearing with Forced Feed Lubrication. The oil 
flow for a bearing with forced feed depends on the number and arrange- 
ment of the oil feed openings and lends itself to analytical c'omputa<- 
tion only in certain special cases one of which is chosen below. 

We assume that the bearing in the plane of symmetry is provided 
with a ring groove in which oil is fed under pressure. We assume 
further, as above, that the pressure along the axis of the bearing is 
expressed approximately by a parabolic type curve the maximum ordinate 
of which is equal to the feed pressure Pq (fig. 16). We shall 
compute the quantity of oil flowing from the clearance in the direction 
of the ends of the bearing, not considering for the present the oil 
flow from the loaded zone. 

According to equation (21) the volume of oil q" flowing per 
second from the two sides of the clearance is expressed as follows: 

„ bh 3 dp 
q . — — • 

6 q dz 

where h is the height of the clearance 
b the width of the clearance 

dp/dz the pressure gradient on the flow direction, that is along 
the Z axis. 

We again turn to figure 3* For the element of width rdqp and height 

h « r^ (1 +l cosqp) 

the flow per second is equal to 

dq" = r ^ 5 dp (33) 

6 q dz 
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By the equation of the parabolic curve 


<lp m 2mp 0 
dz 1 


w 

ro 

vo 


vhere m Is the exponent of the curve, p 0 the feed pressure. Sub- 
stituting this expression in equation (33) and integrating the latter 
between the limits 2n - 0, that is, over the entire periphery of the 
bearing, we obtain the total oil discharge per second: 

«■ . r -^p / (1 ♦ X COB ,)3 d* 

0 

. e, (1 ♦ 1.5 o) - -*> OOo 

3ni 1 qi 

The function (X)* is plotted against X n in figure 17- The latter 
shows that for full eccentricity of the shaft (X = 1) the oil flow 
is two and a half times as large as for central position of the shaft 
(X = 0) . This was to be expected since according to equation (21) 
the flow very much depends on the amount of the clearance. 


The oil discharge per rotation is equal to 

4 ■ - MS <*>0 

H n qnl 

Expressing nr] in terms of (X) (equation (4) ) we obtain 

1“ * (x> 


where 


<3n" = 


13.3 


rVW 0 (X) n 

-H- 2 


(34) 


The function (X) D /(X)* is shown plotted against X in fig- 
ure 18 and with sufficient accuracy is expressed in terms of t]cd/1c ^ 
by the following empirical equation 


1.2 


(X) c 

(X)</inm\ 1,g 
^cEY 2 J 
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Substituting this expression in equation (34), setting m = 2.5 unfl 
expressing (l + d/l)i * 2 by an approximate formula from the binomial 
expansion we obtain 



„ » . 2 r d V ,4 Pog" 
- 

(35) 

where 

•'■-i ( i + i - 2 f) 


The flow is 

»• . 0.4 

(36) 




Seducing expressions (35) and (36) to practically suitable 
we obtain 

units 

l" 

= c^/eec - 2.5 • 10 6 

(37) 

V 

= cm 3 /rotation = 1.6 • 10 ® ^ 50 

kd 0 * 4 A 1 - 2 

( 38 ) 

where A ■ rjn/k 




d the diameter of the shaft in mm 

A the diametral clearance in mm 

2 

Pq the feed pressure in kg/cro 

T) the viscosity in centlpolses 

n the rpm 

k the projected unit bearing load in kg/cm 2 
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1 - the length of the hearing in mm 

<t" = j- (1 + 1.2 d/1) 

Formulas ( 35 ) to -( 36 ) show that the oil flov from the forced 
feed depends to a large degree on the clearance, the ratio l/d, 
is directly proportional to the feed pressure and inversely propor- 
tional to a small power of X . Thus the flow of oil under the 
effect of forced feed, like the flow of oil from the loaded zone, 
is a factor of stability assuring the maintenance of a definite 
value of Tpo/k. 

Range of Application of Formulas ( 35 ) to ( 38 ) Equations ( 35 ) to 
( 38 ) are applicable to the case of forced oil feed from the end of 
the bearing except that in this case the values of q" and q n " 
must be reduced to one fourth. In case the oil is fed through one 
or several holes on the periphery of the bearing equations ( 35 ) to 
( 38 ) may be applied only as a first approximation, .The more holes 
on the periphery and the closer they are arranged in the plane of 
symmetry of the bearing the more accurate the results that may be 
expected from formulas (35) to (38) . In the case of a single hole or a 
small number of holes the actual discharge will be less than that com- 
puted by equations (35) to (38) by an amount which can be determined 
only experimentally. 

If the holes are located in the lower loaded zone of the bear- 
ing (which rarely occurs in practice) there may be expected not 
only a decrease in the amount of oil flowing out but also a change 
of the dependence of the flow on X. The shaft, as Is known, 
changes its position in the bearing with change in X , Its center 
describing a path approximating a semicircle the top point of which 
coincides- for infinitely large value of lpo/k with the center of 
the bearing. This displacement changes the cross-section of the 
oil holes and affects the flow of the oil. 

The grooves and cut-aways in the bearing and shaft as shown in 
fig. 12a fi-nd 12c may increase the flow of the oil more than ten times 
the value computed by equations (35) to (38) and the discharge will 
not depend on X. If the bearing has structural features having 
the object of increasing the oil circulation the flow of the oil 
through the clearance may be neglected and the oil outflow computed 
from the geometric dimensions of the grooves with the aid of the 
usual formulas of the flow of a viscous fluid. 

Total Flow in Bearing with Forced Feed. The total flow of oil 
from a cylindrical bearing with forced feed is made up of the oil 
flowing out of the loaded zone and the oil flowing out under the 
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effect of the feed pressure. As a first approximation we may 
neglect the effect of the oil feed grooves on the pressure 
distribution in the loaded region and consider the total, flow of 
oil from the hearing as equal to the sum of Q* and q" by 
equations (29) to ( 32 ) and ( 35 ) to (38) . The total flow of oil in 
this case is expressed by the equation 

Q cm 3 /sec = q* + q" = 0.8 — — + p 2.5 • lcA-^ ? - Q a (39) 

. 0.6 . 0.6 


where p is a magnitude less than 1. 

Comparison with Test Results. Formulae (29) to (32) and ( 35 ) 
to (3 5 T were compared with the test results on oil flow published by 
Barnard (reference 17 )- The latter determined the flow of oil from a 
full cylindrical bearing of diameter 25.4 mm, 50.8 mm length with 
diametral clearance 0.15 to 0.28 mm. The feed pressure was varied 
between 0.7 to 4.9 kg/cm 2 , the load between 2.9 to 19 kg/cm 2 , the 
speed between 200 to 2000 rpm and the viscosity of the oil between 
12.5 to 43.5 centipoises. The shaft was of hardened steel and the 
bearing of bronze. In the loaded part the bearing had a narrow 
oil distributing groove parallel to the axis of the bearing extend- 
ing up to 5 mm from the ends. The oil flow through the beetring 
was expressed by the "useful pumping coefficient" E equal to the 
ratio of the actual flow of oil per rotation to the volume of the 
bearing clearance : 


E = 2 Qn/itdAL 

The results of the observations were grouped in the form of a set 
of curves of "useful pumping coefficient" as a function of the 
operating parameter rjn/k and the feed pressure Pq (fig. 19) • 
Figure 20 gives the values of the pumping coefficient computed by 
equations ( 32 ) and (38) . Unfortunately Barnard does not show 
separately the values of n, q and k corresponding to the various 
values of X whereas in the right member of our equation there 
enters the value of X in addition to k. This difficulty was 
circumvented by assuming k inversely proportional to X. 

The "structural factor" of equations (32) and (38), equal 
respectively to qn’/X 0, ° and q n "/Xl. 2 # was obtained by dividing 
the values measured by Barnard by X. Hie agreement of figures 19 
And 20 is more than satisfactory. It shows the correctness of 
the fundamental assumptions underlying fo rmula s (29) and (32) and 
( 35 ) and (38) . Although our conclusions regarding the effect of X 
on the oil flow through the bearing are entirely satisfied by the 
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test results of Barnard the effect of the clearance on the oil flow- 
according to our considerations Is considerably greater than accord- 
ing to Barnard who did not at all observe the effect of the clearance 
on the useful pumping coefficient whence the conclusion may be drawn 
that the oil flow according to Barnard is directly proportional to 
the first power of the clearance. Formulas (32) and (38) however 
give a much sharper dependence. Ms point requires special investi- 
gation: 


Comparison of q* and q". We shall consider the relative 
Importance of the flow through the loaded zone q* and that due 
to the effect of the feed pressure q". Dividing equation (28) by 
(32) and dividing out a* and a" which do not greatly differ from 
each other we obtain 
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1.2 


Po 1 
TT k 0.6 


OQ 

ro 
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Substituting the most frequently occurring value of the clearance 
♦ ■ 0.001 we obtain 




x °- 6 


We assume k - 50 kg/cm 2 and shall vary the feed pressure Pq 
between 1 to 6 at. The results of the computations for the various \ 
are given In table 7. Table 7 shows that for a smooth bearing and 
usual relative clearance V » 0.001 the flow under the effect of the 
forced feed considerable exceeds the "natural" flow of the bearing 
especially for small values of \ and large feed pressures. An 
increase in the clearance still further Increases the part played by 
the forced flow. 


Grooves, cut-aways, etc. in the unloaded zone may still more 
sharply change the relation between the "natural" and "forced" flows 
of the oil in favor of the latter and entirely mask the "natural" 
flow. 


Range of Application of Obtained Results. We again mention the 
fact tha t formulas (29) to (32) and (35) to (38) are applicable only 
to the case of smooth shafts and bearings which do not deform under 
load for clearances of the order of f 0.002-0.0005 for constant 
or centrifugal load. In the case of variable or impact 
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load these formulas may he used only as a first approximation which 
is the more accurate the greater the part played by the "forced" 
flow in comparison with the "natural" flow. An analysis of the 
natural, flow of oil for impact and variable loads requires the 
establishment of relations between the eccentricities of the shaft 
and the rate of change of the load, that is extension of the lubri- 
cation theory and its verification by experiment which, incidentally, 
is technically considerably more complicated than in the case of a 
constant load. 


3. HEAT BALANCE OF THE HEARING 

For a steady thermal state the amount of heat generated in 
the bearing is equal to that dissipated. The heat transfer is in 
the following three main parts: (D in the oil flowing out of the 
bearing, (2) in the bush and body of the bearing, f 3) in the shaft and 
parts connected with it. From the body and shaft the heat is removed 
to the surrounding atmosphere by convection and contact. 

The two latter forms of heat transfer lend themselves with 
difficulty to mathematical computation. To a very large extent 
they depend on the special features of the bearing, the surfaces 
of the bearing and shaft, the number and shape of the structural 
details, etc. The heat transfer from the shaft, for example, increa- 
ses sharply if there are rotating masses on the shaft, such as 
pulleys, couplings, propellers, etc. The heat transfer is greatly 
increased on the connecting rod journals of crankshafts in an air 
stream. 

As far as is known to us the only attempt to compute this heat 
transfer (together with the heat given off by the oil) was made by 
Falz (reference 4). T* 10 computation is based on the tests of Lasche 
and is of a very primitive form. Lasche found that the heat transfer 
from, the shaft of a bearing is proportional to the 1.3 power of 
the temperature difference between the bearing and cooling medium. 

The heat transfer equation according to Falz -Lasche has the follow- 
ing form: 


B k cal/hr « adl (t n - t a ) 1,3 


where 

d is the diameter of the shaft 
1 the length of the bearing 
t Q the temperature of the bearing 
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t_ the temperature of the surrounding medium 

0 a coefficient varying with the kind of bearing in the range 

a 50 to 2000 

For bearings of internal combustion engines Falz recommends choos- 
ing the values a » 50 to 2000 depending on the speed of the engine. 

The exponent of (t - t ) which is somewhat greater than 1 showt 
that the most important n part a in the heat transfer of the bearing is 
played by convection and contact and not by radiation which is pro- 
portional to the difference between the fourth powers of the tem- 
peratures . 

In practice a number of cases are known where the heat generated 
in the bearing is given off to the oil and the heat dissipation to 
the surrounding medium may be neglected. To such a case belong for 
example the bearings of an internal combustion engine of the air- 
plane or automobile type enclosed in the crankcase and having a 
high temperature due to the contact with the heated working parts 
of the engine . The heat transfer from the bush of the bearing to 
the crankcase is in this case insignificant due to the anw.11 dif- 
ference in temperature and the heat transfer from the shaft in the 
hollow of the crankcase is likewise small. An exception is pos- 
sibly the case where the crankshaft is directly connected with a 
metal air propeller. 

If we limit ourselves only to the heat dissipation in the oil 
then for a heat balance of the bearing under consideration the 
following equations hold. 

The equation of the heat generated in the bearing: 

R k cal/sec = Pvf =^~ ( — 0.55 . 0 ^ (40) 

427 427 \\J/ k / 

The heat carried off by the oil in unit time: 

R' k cal/sec = c Q — (t r - t n ) (41) 

1000 r 0 


where 

c is the specific heat of the oil in k cal/kg° C 

Q the flow of oil though tjie bearing in cm 3 /sec 
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and t Q "tii® final and initial temperature of the oil. 

c and 7 vary insignificantly with the temperature . The temperature 
of the bearing t n may to a first approximation be considered 

equal to the mean temperature of the oil: 


- 


t^» + t^ 


= t r 


At 
2 ’ 


vhere At is the increment of the temperature in the bearing. In 
this case equation (41) takes the following form: 


E* = 2 c Q 


8 


1000 


(tr 



) 


The relation between the viscosity of the oil and the tempera- 
ture is given by the curve of viscosity against t or by one of 
the many characteristic equations proposed by a number of authors. 
Thus according to Poiseuille (reference 9) : 


1 

A + Bt + Ct 2 


where A, B, C, are constants, characteristic for a given oil. Accord- 
ing to Klesskalt (reference 18 ) : 

„ * -VT 

n = no 15 

where b is a constant . According to Falz (reference 4) : 

11 (O.lt) 2,6 


where i varies from 1.06 for heavy oils to 0.07 for light oils. 
The author found that the viscosity of aviation oils (fig. 21) is 
well expressed by the following relation 


n 


i 

(O.lt) 3 


(42) 


where i = 2.0 to 2. 8 for heavy oils of the type of "brightstock" 
«.nd 1.4 to 1.8 for castor oil and "kastrol" . 


For a steady thermal state the heat given out in the bearing 
is equal to the heat given to the oil: 



32 


NACA 3M No. 1165 


R - E ’> S', P ^ 0.550^ = 2c Q —i (t_ - t Q ) (43) 

427 U k ) iooo n ° 

The simultaneous solution of equation (43) and. the characteristic 
equation of the viscosity permits finding the temperature of the 
hearing for given speed, load and geometric dimensions . 

We shall solve for example the following problem: It is 

required to find the temperature and safety factor of a bearing 
carrying 5000 kg at n = 2000 rpm. The dimensions of the bearing 
are d = 100 mm l/d = 0.8, A = 0.1 mm, feed pressure of oil 3 at, 
temperature of oil 65°, the oil is "kastrol" with i = 1.8. 

The unit load on the bearing is equal to: 

5000 . o 

k = = 63 kg/cnr 

80 


Assume that the temperature of the bearing is 90°. The viscosity 
of the oil is equal to 25 centipoises and A is 


25-2000 

60 


800 


The friction coefficient by equation (18) is equal to 


f = 3.36 • 10~ 9 — ^ — 800 + 0.55 ^ 10 ° • 0.1 = 0.0034 

iooo so 1 - 8 

The heat given off is 

Pvf 5000 * 10.4 . 0.0034 

R = -j-r=-= = 0.43 cal/ sec 

+2T 427 


The oil flow by equation (39) (for P = l) is equal to 

Q = 0.8 


.0.8.2. 2_i, *3.4-o^/ri 


d“ '“4" o ' n + j _ 5 . lo6 , 19 j 0 „3 /sec 


?> 0.6 


d°'V °' 2 


Taking c = 0.5 and y = 0.85 we obtain from equation (43): 

(t„ - 65) 

100000 n 

t n = 91° 


0.43 = 19.7 
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Thus no recomputation is required. The temperature of the oil 
passing out of the bearing is equal to 

65+2 (91 - 65) « 117° C. 

We shall find the factor of safety of the bearing for which 
the author (reference 19) proposes the ratio of the actual operating 
parameter X to the critical parameter X corresponding to the 
start of breakdown of the oil film. The critical value of the 
parameter is found from the limiting thickness of the oil film for 
which there is as yet no contact of the projecting roughnesses of the 
shaft end bearing. This magnitude depends mainly on the accuracy 
of the evaluation and the character of the surfaces of the shaft and 
beaming. We assume that the criticeil thickness of the oil film 
is 0.01 mm which corresponds to the eccentricity: 

0.05 - 0.01 

X = 0.8 

0.05 


From figure 5 for X/d * 0.8 the corresponding value of the magnitude 
T]CD/k\j^ = 0.36 whence the critical X is 

- c 0.36 

X = 9368 . 105 * 335 

10 5 


The safety factor of the bearing is equal to 


x 


-I..*™, 2 . 4 

X kp 335 


The safety factor can be increased in two ways: by an increase 
in the flow of oil through the bearing and by a decrease in the 
temperature of the oil feed. We shall first consider the first 
method. We shall raise the feed pressure to 6 at. The exact com- 
putation gives in this case t « 87°. The viscosity of the oil is 
then equal to 27 centipoises, A increases to 870 and the safety 
factor to 


x 


870 


» 2.6 


335 
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Wo shall now attempt to lover the temperature of the oil 
feed to 50° at Pq a 6 at. In this case the temperature of the 
hearing drops to 84° which corresponds to a viscosity of the oil 
of 31 centipoises, X drops to IlhO and the safety factor to 


335 


M 

ro 


There is noted the insignificant change in the temperature of 
the hearing for relatively large changes in the pressure and the 
temperature of the feed oil. This can be readily explained. By 
feeding a cooler oil ve at the same time increase the friction as 
a result of which the heat generated in the healing increases. On 
the other hand the flow of the oil decreases and the heat dissipa- 
tion drops. A decrease in the inlet oil temperature even to as 
low a value as 20° gives a lowering in the temperature of the bear- 
ing only to 78°. 

A much greater improvement is obtained by increasing the oil 
flow through an increase in the flow cross-sections in the unloaded 
zone, a method which at the same time reduces the friction. We 
assume that a cut-away is made in the hearing of height 0.5 mm on 
the unloaded bearing arc extending over 180°. The friction in this 
region may he approximately computed by the law of Newton: 

T" = t] I ndl = t]U>, f " = T ” = * d ^ 

h 2 2h kid 2h. k 

The friction in the loaded zone of the bearing may be approxi- 
mately computed by the formula of Gumbel preference 5) : 


r = x.7^pi ff 

The total coefficient of friction is equal to 

f « f' + f" - 1.7 + 1 a/^"+ 5^ 

v 1 V k 2h k 

The flow through the hollow (neglecting the effect of the eccentri- 
city of the shaft) is found, by formula (21) : 

rtdh 3 p 

Q = — — 

6r)l 
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The flow through, the loaded zone may he neglected. The temperature 
of the oil is determined from the equation of heat balance (43) . 
Several successive approximations give: 

tn = 65 + 4 = 69° 

The viscosity is then equal to 60 centipoises, A is equal to 
1900 and the safety factor is 

1900 _ _ 

x * = 5.7 

335 


We shall lower the temperature of the inlet oil to 50° . ' The tem- 
perature of the bearing by the equation of heat balance is then 
obtained equal to 58.2°. To this corresponds a viscosity of 97 
centipoises, A is equal to 3100 and the safety factor is 


x 


3100 

335 


9.2 


This example Bhows that by a decrease in the temperature of 
the inlet oil a sharp increase in the safety factor and load bear- 
ing capacity of the bearing may be attained for a single constant 
condition namely an increased oil circulation obtained by structural 
means . If this condition is not satisfied the decrease in the tem- 
perature (and increase in the pressure) of the inlet oil does not 
give any essential result. 

Load Bearing Capacity and Structural Factors. We shall consider 
the effect of structural factors on the temperature and load bear- 
ing capicity of the bearing with forced feed. For simplicity we 
shall assume that the bearing works at a high feed pressure so that 
the "natural" flow of the oil through the bearing may be neglected. ' 
We shall assume that the coefficient of friction varies according 
to equation (19). In this case the heat generated is equal to 

R _ Pvf m Pood 5 /jcA °- 85 

427 2*427 ^0.7' k' 

The oil flow is equal to 

*■ . o.4 

tjTP * 2 

The heat transfer is 


7 

1000 


(*n - *0) 


E* = q" 
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The Increment In the temperature of the bearing Is determined from 
the equation 




/_1A 
V di 


At - const (44) 

i P 0 

The above equation shows that the temperature increment in the bear- 
ing is proportional to the square of the viscosity, the speed and 
the ratio l/d, inversely proportional to the fourth power of the 
relative clearance and ivnersely proportional to the feed pressure . 


We shall now analyze the effect of the heat condition of the 
bearing on its load bearing capacity. We shall make use of the 
expression of Gumbel - Falz (reference 4, 19) based on the Reynolds 
equation for the load bearing capacity: 


k *= const 


1 ) to d^ 
h 0 min Ac 


(45) 


where Iiq . is the limiting allowable thickness of the oil film, 

c « 1 + d/l the correction for the finite length of the bearing. 
Substituting the viscosity formula (40) in equation (45) we obtain: 


k = const 


COd 



(46) 


Substituting in place of At 


its expression from (44) we obtain 


K = const 



(47) 


The above equation shows that the predominating effect on the load 
bearing capacity is that of the bearing temperature. If no special 
means are provided for increasing the circulation of the oil the 
load bearing capacity Is determined principally by the increase in 
the bearing temperature, the analysis of which was given above 
(foimula (44)). 


If forced circulation is used the temperature of the bearing 
is not large and the load bearing capacity is determined mainly by 
the temperature of the inlet oil so that all efforts of the designer 
should be aimed at a reduction of t Q . If the temperature of the 

beetring is given by the individual features of the bearing the load 
bearing capacity only be affected by a change In the structural 
parameters and a choice of the kind of oil used. 


<n 

ro 

VO 
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COMPUTATION CF BEARING 


The computation below refers to the special case where the heat 
given off to the surrounding atmosphere may he neglected and the main 
heat dissipation occurs In the oil. 

The computation reduces to the determination of the tempera- 
ture of the hearing with the aid of the equation of heat balance 
and the characteristic equation of the viscosity. From the tem- 
perature of the hearing is determined the mean viscosity of the oil 
after which hy the known equations of the hydrodynamic theory it is 
determined whether the hearing can sustain the given load and with 
what safety factor. If too small values are obtained for the latter 
an increase in the load hearing capacity is obtained hy one of the 
methods Indicated above. 

The clearance is rarely given as a definite figure hut us ually 
with a tolerance determined hy the manufacturing possibilities. 

This circumstance and also the possible wear of the hearing and 
shaft must he taken Into account In the computation In view of the 
considerable effect of the clearance on the friction and oil flow. 

To Increase the reliability of the computation the oil flow 
and coefficient of friction must he computed for the mini un nn 
tolerance and no wear. The oil feed, however, the cross-sections 
of the oil feed pipes and the output of the oil pump must he 
computed for the marlmim tolerance and with wear taken Into account. 
If the shaft and hearing are of different materials the amount of 
the clearance must he considered In the warm state taking account 
of the difference in the linear expansion coefficients of the 
materials . 

Reference Data. In conclusion we present some data on the 
heat capacities and densities of oils at various temperatures. 

Figure 22 shows the change of the densities of various aviation 
oils with temperature. The density may be determined hy the formula 

7 - 7 0 (1 - a t) 

where 7 q Is the density at the initial temperature, a a constant 
lying between 0 . 0006 * to 0.00070 and on the average equal to 0 . 00068 . 

The heat capacity of the oil Is expressed hy one of the follow- 
ing formulas: 

o _L(o.43 + 0.00081 t), c = 7 = + B (t - 15) 

kg° c V/ V 7 -*-5 
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where 

7 ^ i s the density of the oil for 15° C 

A a constant equal for neutral mineral oils to 0.405 to 0.4-25 
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TABLE I 


Time of publication 

Friction fonnula 

Type of 
lubri - 
cation 

d MM 

1 

a 

t 

Material 

Speed 

range 

m/Bec 

Pressure 

range 

Tempera- 

ture 

°0 

Viscosity 

range 

centipoises 

1885, B. Tower 

f =0.0055 

* 

Oil bath 

101 

1.5 

1/1120 

Steel on bronze 

0.5-2. 4 

7-43 

32 

22-100 

1902, Strlbeck 

0.0016 (l) 2 ^ 3 v l/3 

Ring 

70 

1.95 

1/950 

Steel on babbitt 

1.35-4 

2-22 

20-47 

56-250 

1902, 0. Lasche 

0.0052 ^ vl/ ' 3 
k 

Ring 

260 

0.42 

1/950 

Steel on babbitt 

1-15 

1-15 

28-100 

7-40 

1913, C. Thomas, E. Maurer, L. Kelso 

0.002 ^ 

Ring 

62 

3.95 


Steel on babbitt 



27-43 

40-100 



1-40 

1889, A. Martens 

0.00265 


99 

0.71 

Steel on bronze 

21-41 

30-200 



0.2-0. 5 

1903, A. Kingsbury 

0.00275 (^ff^ 

35 

1.45 

Steel on bronze 

2. 8-2. 4 

32-77 

15-46 



1904, D. Woods and D. Carter 


230 

0.78 

Steel on bronze 

2.3-11 

3.5-10.5 

22-60 

28-120 



1922, D. Stoney, R. posvall, D. Massey 


Bath 

63.5 



Steel on babbitt 

8.5-19 

2. 8-8.4 

38-60 

16-42 

k 0 ' 63 



1922, G. Hovarth and Nelson 
1922, M. Hersey 



50.8 

25.4 

0.23 


Steel on bronze 
Steel on bronze 


28-56 

2.8-17.5 

40-90 

10-30 

38-420 

k 2 / 3 

O.OOH5 









1923, G. Stanton 



25.4 



Steel on bronze 


19 

32-50 

25-209 

kl/4 





1928, D. Goodman 

„ .,0.37 v 0.63 

Bath 

152 



Steel on bronze 

0.75-3.5 

2-15.4 

16-72 

40-54 

Vk 



1930, L. Ulmer 









20-150 












H 


NAGA TM No. 1165 



TABLE II 

Y = 0.001 


Tj 03 

kY 2 

0.01 

0.05 

0.1 

0.5 

1 

5 

10 

50 

100 

X 

10 

50 

100 

500 

1000 

5000 

10,000 

50, 000 

100,000 

It T]Oi 

y * k 

0.00003 

0.000157 

0.000314 

0.00157 

0.00314 

0.0157 

0.0314 

0.157 

0.314 

In o/o of f 

G.3 

25 

39.3 

73.4 

83 

95.7 

97.8 

99.5 

99.8 

0.314 Y 

0.000314 

0.000314 

0. 000314 

0.000314 

0.000314 

0.000314 

0. 000314 

0.000314 

0.000314 

In o/o of f 

63.2 

50 

39.3 

14.6 

8.3 

9.1 

1 

0.2 

0.1 

f" 

0.00015 

0.00016 

0.00017 

0.00026 

0.00034 

0.00039 

0.00039 

0.00041 

0.0004 

In o/o of f 

30.5 

25 

21.3 

12 

8.7 

2.4 

1.2 

0.3 

0.1 


0.00049 

0.00063 

0.0008 

0.0021 

0.0038 

0.0164 

0.032 

0.1577 

0.3147 


TABLE III 

Y = 0.002 


Tj CD 

0.01 

0.05 

0.1 

0.5 

1 

5 

10 

50 

100 

kY 










X 

40 

200 

400 

2000 

4000 

20,000 

40,000 

200,000 

400,000 

It _ TJ(U 

Y ' k 

0.0000628 

0.000314 

0.000628 

0.00314 

0.00628 

0.0314 

0.0628 

0.314 

0.628 

In o/o of f 

6.3 

25 

39.3 

73.4 

83 

95.7 

97.8 

99.5 

99.8 

0.314 Y 

0.000628 

0.000628 

0.000628 

0.000628 

0.000628 

0.000628 

0.000628 

0.000628 

0.000628 

In o/o of f 

63.2 

50 

39.3 

14.6 

8.3 

1.9 

1 

0.2 

0.1 

f" 

0.0003 

0.00032 

0.00034 

0.00052 

0.00068 

0.00078 

0.00078 

0.0008 

0.0008 

In o/o of f 

30.2 

25 

21.3 

12 

87 

2.4 

1.2 

0.3 

0.1 

f 

0.00099 

0.00126 

0.00150 

0.00429 

0.00759 

0.0328 

0.0642 

0.315 

0.629 
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TABLE IV 


\|f = 0.0005 


3 m 

0.01 

0.05 

0.1 

0.5 

1 

5 

10 

50 

100 

k\|/? 










X 

— 

— 

25 

125 

250 

1250 

2500 

12,500 

25,000 

it _ i](D 

y ’ k 

— 

— 

0.000157 

0.000785 

0.00157 

0.00785 

0.0157 

0.0785 

0.157 

In o/o of f 

— 

— 

39.3 

73.4 

83 

95.7 

97.8 

99.5 

99.8 

0.314 \|r' 

— 

— 

0.000157 

0.000157 

0.000157 

0.000157 

0.000157 

0.000157 

0.000157 

In o/o of f 

— 

— 

39.3 

14.6 

8.3 

1.9 

1 

0.2 

0.1 

f" 

— 

— 

0.000085 

0.00013 

0.00017 

0.00019 

0.0002 

0.0002 

0.0002 

In o/o of f 

— 

— 

21.4 

12 

8.7 

2.4 

1.2 

0.3 

0.1 

f 

— 

— 

0.0004 

0.00107 

0.0019 

0.0082 

0.016 

0.788 

0.1573 


TABLE V 

V = 0.001 


Tim 

0.01 

0.1 

1 


1000 

ck y 






Const Tjy 2 

0.0000314 

0.000314 

0.00314 

0.0314 

0.314 

o/o of B 

5.4 

36 

85 

98.39 

99.83 

Const kv 

0.00055 

0.00055 

0.00055 

0.00055 

0.00055 

In o/o of E 

94.6 

64 

15 

1.7 

0.17 

B 

0.0005814 

0.00096 

0.00369 

0.03195 

0.3145 
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TABLE VI 

= 0.002 


n g> 

2 

LT 

0.01 

0.1 

i 

10 

100 

Const rjv^ 

0.0000157 

0.000157 

0.00157 

0.0157 

0.157 

In o/o of E 

1 . 14 

12.5 

59 

93.5 

99.3 

Const kv 

0.0011 

0.0011 

0.0011 

0.0011 

0.0011 

In o/o of E 

93. 6 

87 . 5 

41 

6.8 

0.7 

E 

0.001115 

0.001257 

0.00267 

0.0168 

0.1581 


TABLE VII 



X 

250 

500 

1000 

5000 

10,000 


f = 1 at - 

5.5 

3.3 

2.55 

0,95 

0.6 

q M 

at. 

17 

11.5 

7.7 

2.9 

1.8 

v 1 

L 5 at . 

27.5 

19 

12.4 

4.8 

3 
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Fig. 10 
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